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Abstract 

In this paper we provide three contributions to the field of channel sounding waveform design in asynchronous 
Multi-user (MU) MIMO systems. The first contribution is a derivation of the asynchronous MU-MIMO model and 
the conditions that the sounding waveform must meet to independently resolve all of the spatial channel responses. 
Next we propose a chirp waveform that meets the constraints and we show that the MSE of our system meets 
the Cramer-Rao Bound (CRB) when the time offset is an integer multiple of the sampling interval. Finally we 
demonstrate that the channel capacity region of the asynchronous system and synchronous system is equivalent 
under certain conditions. Simulation results are provided to illustrate the findings. 

Index Terms 

Multi-user, MIMO, Asynchronous, Channel sounding, chirp waveform 

I. Introduction 

MULTIPLE input multiple output (MIMO) communication is a powerful way to increase the channel 
capacity and combat the effects of fading in wireless communication systems [1], [2]. The channel 
state information (CSI) is very crucial for data detection and channel equalization in MIMO systems, 
as well as for determining the capacity of a given MIMO deployment configuration. Multi-user MIMO 
(MU-MIMO) is being deployed in both the LTE advanced and WiMax standards [3], [4] as a method for 
spatially multiplexing multiple user data streams in the same time/frequency slot in order to maximize 
channel utilization. 

In order to ascertain the capacity and hence the effectiveness of MU-MIMO, it is crucial to perform 
channel sounding tests to determine the spatial channel response across all of the MIMO channels. Channel 
sounding is much more compUcated in MU-MIMO systems when trying to sound the channel from 
multiple sources that have significant distance separation. Furthermore, one or all of the users may be 
moving relative to each other. 

Most of the MIMO channel sounders [5]-[8] simply extend the single input single output (SISO) channel 
sounders, which use a multiplexing switch to scan the different antennas sequentially. The sequential 
techniques avoid the cost and complexity of the multiple receiver architecture. The works in [9]-[12] 
proposed a semi- sequential architecture by using parallel receivers; however, the switching is still expected 
at the transmitter side. Any system that uses switching will be unable to capture true cross channel statistics 
estimates, which are required in order to estimate channel capacity and other relevant channel metrics. 
This shortcoming is usually ignored, or explained away by stating that the channel coherence time is 
longer than the time it takes to cycle through all channels, but as the size of MIMO arrays get larger, 
such coherence-time constraints become increasingly unrealistic. Another drawback to systems that have 
to cycle through antennas is that they are limited in the Doppler resolution that can be resolved. Doppler 
resolution is inversely proportional to the channel sounding waveform repetition rate. In switching systems, 
this repetition rate increases with the number of antennas inherently limiting the Doppler bandwidth that 
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Fig. 1. Multi-user MIMO systems 



can be resolved. In this paper, we propose a code-base sounding waveform that enables both increased 
Doppler resolving power as well as the ability to estimate channel cross correlation statistics. 

Ideally, each of the transmitters and receivers would be time- synchronized to the same atomic clock 
timer. Under certain test scenarios, this may be feasible and even straightforward to accomplish. How- 
ever, given that channel sounding is usually performed with signal generators at the transmitter and 
signal analyzers at the receiver, there is no easy duplex channel to facilitate clock synchronization. GPS 
synchronization may be used, but is not available for indoor channel sounding. Of course, even when it 
is possible to perform synchronization across all nodes, this will come at the cost of increased expense 
and complexity. 

This paper addresses waveform design for channel sounding in such an asynchronous MU-MIMO 
system. We provide the following contributions: i) We derive the asynchronous MU-MIMO model and 

we show the conditions that the sounding waveform must meet in order to independently resolve all of 
the spatial channel responses, ii) We propose a chirp waveform that meets the constraints. We find the 
Cramer-Rao Bound (CRB) for any channel sounding system and then show that the MSE of our system 
meets the CRB when the time offset is an integer multiple of the sampling interval. For fractional time 
offset, the performance is degraded by an unknown pulse shaping matrix; iii) We demonstrate that the 
capacity is invariant to the channel collected via asynchronous collection versus perfectly synchronous 
collection as long as either the transmitter side or the receiver side is synchronized. 

II. System Model 

We consider an MU-MIMO system with Mt transmit nodes and receive nodes. Assume Nt and 

A^,. to be the total number of antennas for all the transmit nodes and all the receive nodes respectively, 
then the system has x A^,. channels to be sounded, as illustrated in Fig. 1. Each channel is assumed 
to be frequency-selective and block-wise time invariant over the duration of the sounding waveform. We 
denote the channel response by /ii,m[^] for the multipath channel response between the ith transmit antenna 
and the mth receive antenna at lag /, where i = 1, 2, . . . , iV^, m = 1,2,..., iV,., and Z = 0, 1, . . . , L — 1. 
The continuous clock mismatch delay of the link between the i*'' transmit antenna and the m*'* receive 
antenna is denoted by Ci,m- In the asynchronous multi-user system, each transmitter and receiver has one 
local oscillator (LO) which is not synchronized between nodes; thus the clock mismatch delays will be 
equivalent if transmitter and receiver antennas belong to the same transmitter-receiver pair but are different 
otherwise. Let T denote the sampling interval. By dividing the clock mismatch delay Ci,m by T, we can 
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Fig. 2. Illustrations of Multi-user MIMO channels with synchronous and asynchronous transmissions 



obtain 

where di^^ denotes the integer part and jii^rn £ denotes the fractional part. Without loss of 

generality, we let /ij ,„ be positive in this paper; i.e., yUj „ ^ (0, \] - In the following, we will use parentheses 
(•) to represent the continuous-time function and use square brackets [•] to represent discrete-time function. 



A. With only integer offset 

We first consider the scenario with only integer offset. We model the effects of integer offset as part of 
the channel response. The first di^m taps of the channel response are considered to be zeros. Without loss 

of generality, we assume that rfi „ = 0, m = 1, . . . , A^^- Assume that the number of nonzero taps of the 
channel response between the i*'* transmit antenna and the m*^ receive antenna is Lj ,„ and Lj.m+'^j,m < L. 
Two examples of MU-MIMO channels with synchronous and asynchronous transmissions, respectively, 
are illustrated in Fig. 2, where Mt = 3, A^* = 3, Mr = 1, A^^ = 1- 

Each transmit antenna of the channel sounder sends out a signal that consists of periodically repeated 
sequences [13]. Let Si be the cyclic sounding waveform with period N transmitted by the i*'^ antenna. 
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Fig. 3. Illustration of channel sounding in Multi-user MIMO systems 



The signal at the receive antenna m at time n can be expressed as 

Nt L-1 

i^m[n\ = XlXl hi,m[n\si[n - k] + Zm[n\, n^O,...,N -1, 



=1 k=0 



(2) 



where Zm[n] is additive noise at the receive antenna m at time n assumed to be white Gaussian (AWGN). 
We can rewrite (2) in matrix form as 



Nt 
1=1 



(3) 



where and have dimension Nxl, hj^^ is of dimension L x 1, and Sj is ToepUtz matrix of dimension 
N X L defined as 



Si[0] 

Si[L - 1] 
. Si[N-l] 



Si[l - L] 
Si[2 - L] 

si[0] 
Si[N - L] 



(4) 



Matched filters are used at the receivers to estimate the channel response, as illustrated in Fig. 3. For 
estimating the channel response hi^m, the matched filter in matrix form is Sf . We can express the channel 
estimation as 

^i,m ~ ^f^m — ^f^i^i,m + S^S^h„^^ + Sfz^ (5) 
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The SfSi term in (5) is the auto-correlation matrix for the i*'* channel sounding sequence 



Sf^Sj 
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Ri[L - 2] 



(6) 



_Ri[-L + l] ■■■ Ri[0] 
where -Rifr] is the periodic autocorrelation function for the i*'* waveform, which is defined as [14] 

N-T-l N-1 

Ri[r]^ Yl sMs*dn + r]+ J] sMs> + r-N]. (7) 



n=0 



n=N-T 



The Sf^Sy term in (5) is the cross-correlation matrix between the i*'* and the v*'* channel sounding 
waveforms, which can be written as 



Civ\L 2] 



(8) 



where Cro[r] is the periodic cross-correlation function between the itu and Vth sequences, which is defined 
as [14] 



N-T-l 



N-1 



Civ[r]= Yl Si[n\sl[n + t] + ^ Si[n]sl[n + r - N]. 



(9) 



n=0 



n=N-T 



B. With fractional offset 

Next, we consider the scenario with fractional offset. Due to the existence of a fractional offset, pulse 
shaping filter has to be taken into account here. Let g{t) denote the pulse shaping filter defined on the 
interval [— MT, MT], the transmitted continuous-time waveform at the i*'' antenna Si{t) can be expressed 
as 



^ii^)= Y Si[n]g{t - nT). 



(10) 



At the receiver, the n*" sample at the receive antenna m can be expressed as 

^mi^] = r^X^T + Hi^rnT) (11) 

Nt L-1 

= [^] Si{nT + fXi,mT - IT) + Zm [n] 

1=1 1=0 

Nt L-1 oo 

= I] 5] ^^.-[^] Yl ^i\.k\9(nT + iJii^rnT -IT- kT) + z^H, n = 0, . . . , TV - 1. 

1=1 Z=0 k=—oo 

Let y — n — k. Since g{t) is only defined on the interval [—MT, MT\ and jii^rn £ (0, \\, we can obtain 

-MT <yT + fXi,^T-lT <MT. (12) 
Since < / < L — 1, we can further obtain 

-M<y<M + L-2. (13) 
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Plugging y — n — k into Eq. (11), we can rewrite (11) as 

Aft L-1 M+L-2 

r^H = Y.ll ^^[^ - y\9iyT + '"^.-^ - lT)h,,m[l] + Zm[n], n = 0, . . . , TV - 1. (14) 

i=l Z=0 y=-M 



The Eq. (14) can be rewritten in matriK form as 



Aft 



(15) 



where = [r^[0], r^[l], • • • ,r^[N — 1]]^ has dimension x 1, Sf is Toeplitz matrix of dimension 
N X (2M + L-1) defined as 



Si[M] 
Si[M + l] 



Si[-M - L + 2] 
Si[-M-L + 3] 



(16) 



_Si[M + N-l] ■■■ Si[N-l] ■■■ Si[-M -L + N + l] _ 
and G{iJ,i^rn) is Toeplitz matrix of dimension (2M + L — 1) x L defined as 

^(-MT + /i,,„T) ••• ^((-M-L + l)T + /i,,^T) 



gi{-LT + Ili^rnT) 

gi{-{L - 1)T + fii,mT) 
gi{-{L-2)T + f,i,^T) 



(17) 



5((M + L-2)T + /Xi,^T) ■■■ gi{{M - 1)T + i^i^mT) 
By left multiplying received vector with the matched filter (Sf the output can be expressed as 

Km = (Sf )^r4 = (Sf )^Sf G(/.,^)h,„ + X;(Sf )^Sf G(/., Jh,,„ + (Sf )^z^, (18) 

where (Sf )^Sf and (Sf )^Sf are of dimension (2M + L - 1) x (2M + L - 1) expressed as 



(19) 
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III. Objective of Sounding Waveforms Design 

In this section, we show the conditions that the sounding waveform must meet in order to independently 

resolve all of the spatial channel responses. 

Firstly, with only integer offset, the estimated channel response in (5) consists of both the desired signal 
term Sf Sjhj^m and the error term J2v^i ^f^viiv,m + Sf z. The ideal situation would be Sf Sjhj^^ = hj^^ 
and Sf S^h^ „i — 0. Hence, the objective is to make the waveforms satisfy 

SfS,= ('-- ' = " (21) 
which is equivalent to the condition that 

10, if < \t\ < L 

Ciy[T] = 0, if -L<T<L (23) 

It is worthwhile to point out that the criteria for an asynchronous MU-MIMO system in (22) and (23) are 
stricter than those for a synchronous MU-MIMO system. In the synchronous MU-MIMO system, since all 
the transmitters and receivers are synchronized, it is unnecessary to model the effects of clock mismatch 

delays as part of the channel response. Hence, the overall length of chaimel response L can be set as 
small as L = max{ Lj„i}. However, in the asynchronous MU-MIMO system, we have to take the clock 

i,m ' 

mismatch delays into account, and the minimum L can be set as L = max{Ljm + di m}. 

i,m ' ' 

Secondly, with fractional offset, the interference from other channels ^t,^j(S,f )^Sf G(/Lii.m)h.u m in 
(18) can be eliminated by making (Sf)^Sf = 0(2m+l-i)x(2M+l-i)- i ^ v. However, h^m can not be 
obtained from the signal term (Sf )^Sf G(/Xi_m)hj^^ since \ii^rn is unknown. Let the waveform satisfy 

^gF^i/gF ^ ) I{2M+L-l)x(2M+L-l), if « = ^^4) 
* " (^0(2M+L-l)x(2M+L-l), if %i^V 

which is equivalent to the condition that 

i?4r] = r' , (25) 

^-^ |t|<2M + L-1, andr^O 

if |r|<2M + L-l (26) 

For comparison purposes, we list various scenarios and the constraints they impose on the optimal 
sounding waveforms in Table I. 

Additionally, since a high peak-to-average power ratio (PAPR) of a transmitted waveform dramatically 
reduces the efficiency of the communication system, we desire low PAPR waveforms. The PAPR of the 
waveform Si is defined as 




max I Si [n] 



|2 



PAPR(si) = . (27) 
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TABLE I 

Constraints on optimal sounding waveforms for various scenarios 



Configuration 


Auto correlation requirements 


Cross correlation requirements 


SU-SISO [15] 




N/A 


SU-MIMO or 

Synchronous MU-SISO or 
Synchronous MU-MIMO [16] 


(1, r = 
= { 0, |r| < max{Li,^} 


Civlr] = 0, T < max{Li,^} 


Asynchronous MU-SISO or 
Asynchronous MU-MIMO 
with only integer offset 


Jl, r = 
^'^'^^^ = ] 0, |t| < max{Li,m + di,m}, 

\ i,7n 


Ci„[T] = 0, |t| < max{Li,m 4- di,m} 

i,m 


Asynchronous MU-SISO or 
Asynchronous MU-MIMO 
with fractional offset 


jl, r = 


Civ[T] = 0, |r| < max{Li,^ -|- di,^} + 2M-1 



IV. Sounding Waveform Design 

To achieve minimum mean square error (MSE) in the channel sounding for asynchronous MU-MIMO 
systems, the waveforms transmitted from the multiple antennas must have impulse-like auto correlation 
function and zero cross correlation as indicated in (22), (23), (25) and (26). The zero cross correlation 
requirement is due to the possibility that the mismatch delay, rfj may be very large, i.e. the maximum 
of Lj m + di^m is essentially unbounded. In [16], [17] a constant-magnitude polyphase channel sounding 
sequence was proposed which has minimum PAPR but does not satisfy the minimum MSE condition. 
The advantages of using chirp waveforms, including the low PAPR property, have been recognized in 
semi-sequential MIMO sounders [9], [11], [12] and other channel sounders [18]-[24]. In this paper, we 
construct the chirp waveform set to satisfy the minimum MSE conditions for asynchronous MU-MIMO 
channel sounding. We design the channel sounding chirp waveforms as 



[n] = exp { j27r^(n + + 2)} , n = 0, . . . , TV - 1, (28) 



where j = \/^. We restrict both the p and N to be powers of 2 and N > 2p, thus, the ratio N/ (2p) is 
always even. The auto-correlation R'^^t] and cross-correlation C^^'^^fr] are proved in the appendix to be 

{1, if r = wf 
-1, if T = {2w + l)g (29) 
0, otherwise 

w^O,...,{p-l) (30) 

C^P'1)[t] = 0,t = 0,1,...,N -1 (31) 

To satisfy (22), we require 

N > 2p^^L- (32) 

To satisfy (25), we require 

iV>2p,„ax(2M + L-l), (33) 

where pmax is the largest value of p for the set of MU-MIMO channel sounding waveforms. 
With only integer offset, the estimation of channel response in (5) becomes 

= hi,m + Sf z„ (34) 
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The mean square error at the receive antenna m can be obtained as 



(35) 



We assume that AWGN z„j has auto-correlation matrix R^,^ = £'{zmZ^} = 2(T^Ijv. Hence the noise term 
Sf z has auto-correlation matrix = E{Sfzz^Si} = 2(7^1^. The MSE^ in (35) can be further derived 

as 



MSE„ 



\\S^z„,f \ = El z^S.Sf z^ ^ = 2La^^ 



We have derived the Cramer-Rao bound (see Appendix) as 



(36) 
(37) 



The results from (36) shows that the MSE of our system meets the CRB. 
With fractional offset, the output of the matched filter in (18) becomes 

hj„, = G(^,,„)h,,„ + (Sf )^z,„. (38) 

We can see that the channel response vector hi^m is multiplied by the pulse shaping matrix G{iJ,i^rn), 
which is a function of the fractional offset iii^m- The noise term (Sf )^z^ has auto-correlation matrix 
R^ = £;{(Sf )^z„z^Sf } = 20-2I2M+L-1. Given the observations hf^, which is a (2M + L - l)-length 
vector, we need to jointly estimate /Xj^^ and L-length channel response vector hj „j. The estimation of 
and Yii^rn is given by 

{/I, hi,„,} = arg min| |hf^ - G(/x)h| I2. (39) 

/i,h ' 

The iterative solution of Eq. (39) can be expressed as 



.(n+l) 



— are mini I hf! 



(n+l) 



1 



G 



H 



(n+l) 



(40) 



(41) 



Since g{t) is a differentiable function, Newton's method can be employed to solve the minimization 
problem in Eq. (40). 

A. Channel Estimate Averaging 

The above analysis only considers at most the (2M-I-L — 1) -length matched filter output. As is apparent 
from (29), each matched filter output contains 2p copies of the channel impulse response estimate. Based 
on this, a natural question arises: can the MSE of the channel estimation be decreased by averaging these 
2p outputs? Next we will show that averaging the multiple segments of the matched filter output can not 
refine the above results and each branch has the same channel sounding performance. It is possible to 
drive down the MSE by averaging the outputs from multiple sequences as long as the channel is static 
across the soundings. We only show that inter-symbol averaging can not decrease the MSE. 

If we assume the channel is static for at least A^^-length matched filter output, we can write the A^^-length 
output as 



i,m 



Si[M+l] 
Si[M + N -I] 

Ri[-i\ ■■ 
Ri[-N + l\ ■■ 



Si[M+l 
Si[M + 2 



N] 
N] 



H 



(42) 



Si[M] 

Ri[2M + L-2] 
R-[2M + L-Z] 



Ri[2M + L-N -l\ 



G(/Xi,„)hi 
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where denotes the noise part of the A^-length output. 



Si[M + l] 



Si[M + 1-N] 
s-[M + 2-N] 



-, H 



(43) 



_s,[M + N-l] ■■■ Si[M] 
By substituting (29) into (42), the A^^-length output can be further written as 



h , 



••2p ' ' 



(44) 



We can see from (44) that the iV-Length output hf^ consists of 2p repHcas of the channel response 
G{iJ,i^m)^i,m plus the fihered noise z^. If the noise z^ were white, we could increase the SNR by a 
factor of 2p by averaging the 2p segments of the matched filter output. However, the noise z^ has the 
auto-correlation matrix 



(45) 



= 2al 



Ri(-i) 

Ri{-N + 1) 



Ri{N - 1) 
Ri(N - 2) 

R-{0) 



which is not a multiple of the identity matrix. In addition, according to (29) and (45), the noise z^ 
has auto-correlation function values R[N/2p] = —2a'^ and R[N/p] = 2af^^. Therefore, averaging the 2p 
output segments with the distance N/2p can not increase the SNR or improve the MSE in (36). The above 
analysis indicates that the performance of channel sounding is independent of the parameter p adopted 
by the waveform. 



B. Capacity Analysis 

Let hi^mif) denote the channel frequency response at frequency / of the synchronous MU-MIMO 
system and hi^mif) denote the channel frequency response at frequency / of the asynchronous MU- 
MIMO system. The equations hi^rn{f) and him(f) should satisfy 



(46) 



The matrix form of the channel frequency response at frequency / for the synchronous MU-MIMO and 
asynchronous MU-MIMO can be written as 



H(/) 



H(/) = 



Kiif) hiif) 

IhiMf) ■■■ 

hiAf) kiif) 

hiM) ••• 



hNt,i{f) 

hNt,Nr if) 



(47) 



h 



Nt,N, 



if) 



(48) 
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The capacity of the synchronous MU-MIMO channel and synchronous MU-MIMO channel are given by 

1 



c — 

^syn — ^ 



Casyn — ^ 



log2det(l+^H^(/)H(/))d/ 
log2det(l+|.H^(/)H(/))d/, 



(49) 
(50) 



where p denotes the signal-to-noise ratio (SNR) and B is the signal bandwidth. 

If the Mt transmitters are synchronized or they share a single LO, the clock mismatch delays should 
satisfy 

f^m = Cl,m = Cz.m = • • • = C,Nt,rm Tfl = 1,2, ... , Nr. 



(51) 



From (46), we can have 

H(/) 
Thus, we can obtain 



e-^^^'f"^ 




H(/) ^ AH(/). 



(52) 



H^(/)H(/) = H^(/)A^AH(/) 
= H^(/)H(/). 



(53) 



If the Mr receivers are synchronized or they share a single LO, the clock mismatch delays should 
satisfy 

l-i — Ci,l — Ci,2 — ■ ■ ■ — Ci,Nr, ^ = 1, 2, . . . , A^t- 



(54) 



From (46), we have 



H(/) = H(/) 



g-j27r/.i 

e 



Q p-32-nfb2 







3-j27r/tAr 



= H(/)B. 



(55) 



Thus, we obtain 

H^(/)H(/) = B^H^(/)H(/)B. 
Since B is a unitary matrix, we have 



det(I+^H^(/)H(/)) = det(l + ^B^H^(/)H(/)B 



P -dH-TtHi 



= det(l + |.H^(/)H(/)). 



(56) 



(57) 



From (49) and (50), we can conclude that the capacity region is the same as long as one side has a single 
LO. However, if both sides have multiple LOs, 



det ( I + ^H^(/)H(/) ) ^ det I I + ^H^ (/)H(/) 1 , 



(58) 



the capacity region is not the same. We need to synchronize the LOs at either the transmitter side or 
receiver side to determine the capacity. 
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Fig. 4. Multipath channels used in the simulations 



V. Simulation 

We set up the simulation scenario as follows: Mt = 2, Nt = 3, Mr = 3, Nr = 3. Thus, we have 9 
channels to be generated. We set the integer offset for the first user to be zero and the second user to 
be 5. All the 9 channels have 10 nonzero taps. The AWGN is assumed to have the same SNR = 25 
dB at the 3 receive antennas in the simulation. To make sure that the sounding waveform satisfies the 
criteria in (33), we use the parameters as: N = 128, p = 1, 2,4 for the three waveforms respectively. The 
generated channel responses are shown in Fig 4. Fig 5 shows the magnitude of the A^-length outputs after 
the matched filter. We can see that the iV-length outputs include 2p = 2,4,8 segments of the sounded 
channel responses. The MSE of channel sounding (generated by 10000 Monte-Carlo runs) and CRB are 
presented in Fig 6. The result with only integer offset matched CRB very well. The MSE of the proposed 
waveform with fractional offset is plotted as well. We randomly generated the fractional offset in each run. 
The performance is degraded by the unknown pulse shaping matrix. The MSE of the waveform presented 
in [16] is also simulated (10000 Monte-Carlo runs) and compared with the proposed waveform. We can 
see that the proposed chirp waveform offers a significant MSE improvement. 

VL Conclusions 

In this paper, we derived the conditions that the sounding waveform must meet in order to independently 
resolve all of the spatial channel responses in the MU-MIMO system. We identified a chirp waveform 
that meets the constraints, then we showed that the MSE of our system can meet the CRB. Finally, we 
demonstrated that the capacity region of the asynchronous system is the same as the perfectly synchronous 
collection as long as either the transmitter side or the receiver side is synchronized. This allows for a 
channel sounding system where only one-sided synchronization is necessary allowing for lower-cost and 
lower-complexity MU-MIMO channel sounding. 
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Fig. 6. MSE of channel sounding 
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Appendix A 
Proof of correlation of chirp waveform 

1. Auto Correlation R^^[t] 

N-T-l 

i?(p)[r] = J2 s'-P\n]s*^^[n + T] (59) 

n=0 



where 



n=N-T 
N-l 

n=0 

n=0 ^ 

exp I - j2TT^{n + 1 + T)(n + 2 + r) 
= — ^exp jj27r-^(-2rn-r^-3r) I 

n=0 ^ 

= ;^exp|j27r-^(-T2 -3r)| ^ (^expj -j2 



^ ^ 2pT\Y (N, T^k"" 



V exp -J27r^ ' "'^^ (60) 

^ V I N \ ) |0, otherwise 

fc = 0,...,(2p-l) (61) 

If r = kN/{2p), i?(p)[r] = exp{jA;7r(-A;iV/(2p) - 3)}. Because kN/{2p) is even, i?^^') [r] = 1 if A; is even 
and i?(^)[r] = —1 if A; is odd. Therefore, 

{1, if r = u)^ 
-1, if T = (2w + l)g (62) 
0, otherwise 

= 0, . . . , (p - 1) 
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2. Cross Correlation C(P'«)[t] 

N-T-l 

C(P'1\t] = Yl s^''\ri])s<'i\n + t] + ^ s^P\n]s*^'i^[n + t - N] (63) 



n=0 n=N—T 
N-l 



n=0 

n=0 ^ ^ 



expj -j27r-|:(n+ l + r)(n + 2 + r) 

n=0 ^ 

— (3p — 3g — 2qT)n + (2p — 2q — Sqr — qr'^)] 



\CM[t]\'^ = C(P''?)[r]C(^'^)*[T] (64) 
1 ^-'^-^ ( 1 

n=0 m=0 ^ 

— (3p — 3g — 2gr)(n — m)] 

iV-l iV-l 

N 



n=0 Tn=0 ^ 



exp <j j2n^[{p -q){n + m) 



-(3p-3?-2gT)] 
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Let 5 — n — m 



N-lN-l . . 

n—n x—n V ' 



where 



Thus, 



n=0 5=0 



exp { j27r— [(p - q){2n - 5) 



-{3p-3q-2qT)] 

1^-1 r 1 

<5=0 

+ {3p - 3q - 2qT)6] 

n=0 ^ 

+(3p-3g-25T)5]| 



iV-1 



n=0 



^ ( exp <^ j27r 



2{ p-q)S 
N 



0, 



5 = 

2q 

Otherwise 



u = Q,...,{2q-l) 



2q-l 



\C^'''W = ^Eexp|j27rl 



M=0 



- (P - 



+ (3p — 3g — 2qT)u 
-}_^exp<^j27r 

M=0 

+j7rM 3(- — 1) 



N 
2q 



u 



1)- UT 

'Aq 



(65) 
(66) 

(67) 



Without the loss the generaUty, we let p > q. Because p and q are both the powers of 2, p/q is an 
even integer. Because N > 2q and N is the power of 2, N/{Aq) is an integer. So the term [—u^{p/q — 
l)N/{Aq) — ut] is an integer as well. Thus, 



^ 2g-l 

- ^ exp <; jTTM 



AT 



u=0 



3(^-1) 



(68) 
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Because p/q is even, the term [3{p/q — 1)] is always odd. Thus, 



2q-l 



(69) 



u=0 



Thus, 



C(P''?)[r] = 0, r = 0,l,...,Ar-l 



(70) 



Appendix B 
The Cramer-Rao Bound 

The in (3) is Gaussian distributed as 

- 7V^/x(hi,^),2a^I^^ 

where //(hi^m) = Sjhj^^ + Ylj^^i ^j^j,m- The Fisher information matrix can be obtained as 



dh,. 



sfs, 



2< 



(2<Iiv)- 



(71) 



(72) 



So the Cramer-Rao Bound is 



2al 



var(hi,„) > 2La^ 



(73) 
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